Thermal noises and noise compensation in high-reflection multilayer coating 
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Thermal fluctuations of different origin in the substrate and in the coating of optical mirrors 
produce phase noise in the reflected wave. This noise determines the ultimate stabilization capability 
of high-Q cavities used as a reference system. In particular this noise is significant in interferometric 
laser gravitational wave antennas. It is shown that simple alteration of a mirror multilayer coating 
may provide suppression of phase noise produced by thermorefractive, thermoelastic, photothermal 
and thermoradiation induced fluctuations in the coating. 



I. INTRODUCTION 

In different optical systems phase noise and shot noise 
are usually considered as the main fundamental limit- 
ing factors of sensitivity. These noises determine the 
standard quantum limit of classical gravitational wave 
antennas [1]. At the same time in laser systems cou- 
pled with high-finesse optical resonators fundamental fre- 
quency stability may be determined also by other fun- 
damental effects connected with mechanical, thermody- 
namical and quantum [2] properties of solid boundaries. 
Many of these effects firstly identified and calculated on 
the forefront of laser gravitational wave antennas [3-8] 
are becoming increasingly important in other high-Q op- 
tical systems [9-13]. 



II. A NOISY ZOO 

Excess optical phase noise is added to the wave re- 
fiected from the mirror of the optical cavity due to varia- 
tion of boundary conditions produced by fluctuations of 
the surface and refractive index in the surface layer. 

I. Substrate Brownian noise. Historically the first 
noise of this origin identified as a problem in gravitational 
wave antennas was thermal intrinsic noise produced by 
internal friction which may be related to the surface us- 
ing fluctuation-dissipation theorem [3] . This noise is fre- 
quently noted as "Brownian motion" or simply thermal 
noise. 

The spectral density of surface fluctuations of a mirror 
averaged over the Gaussian beam radius w (determined 
at 1/e^ decay of intensity) is the following: 
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where fc^ is the Boltzmann constant, T is the temper- 
ature of the mirror, 0(/) is the loss angle of the mirror 
substrate at a frequency / of analysis, as is the Poisson 
ratio and Yg is the Young's modulus of the substrate. 

II. Coating Brovifnian noise. In the same way in- 
ternal friction in the material of the coating produces ad- 
ditional fluctuations, which may be significantly stronger 
due to high level of acoustical losses in the material of 
the layers even though the thickness of this coating is 



relatively small. Effective thermal noise associated with 
the coating is determined by the same relation (1), with 
[8, 14]: 
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Here and below averaging over the multilayer coating is 
performed as follows: 



Xidi + XiJh 
di + dh 



(3) 



Indexes I and h denote constants for low and high re- 
flection coating layers with thicknesses di and c?/i, and 
djsi = N{di + dh) is the total thickness of the coating 
with N double layers. This noise has been experimen- 
tally measured in specialy designed Thermal Noise In- 
terferometers (TNI) and may limit the sensitivity of the 
next generation of gravitational wave antennas [14]. 

III. Substrate thermoelastic noise. In [4] Bragin- 
sky, Gorodetsky and Vyatchanin (BGV) suggested that 
temperature fluctuations in the bulk of the mirror trans- 
formed through thermal expansion to surface fluctuations 
should produce additional mechanical noise. Using both 
the FDT approach and the Langevin method for the anal- 
ysis of thermal correlation functions analogous to that 
developed earlier by van Vliet [15, 16] they have shown 
that this noise may be understood as the noise produced 
by thermoelastic damping. 
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where as is the coefficient of thermal expansion, Cs is 
specific heat capacity, Ks is thermal conductivity, and ps 
is the density of the substrate. 
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IV. Coating thermoelastic noise. The same as III 
but noise is caused by thermal expansion of the coating 
or from another point of view by coating thermoelastic 
losses [6, 7]: 



Qcoat _ ^ksT'^jl + as)^ald% -.coat /r\ 
Otc = 

2 sinh ^ — sin ^ + i?(cosh ^ — cos ^ ) 

'B^ cosh i + cos ^ + 2R sinh ^ + i?^ (cosh ^ - cos ^) 



^=^2cod%{pC){K-^) R. 



(pC) 



For silica-tantala coating ac is within 10% from sim- 
pler (a) averaged as (3) and for ^ <C 1 the following 
approximation is valid: 
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V. Coating thermorefr active noise. In [5] BGV 

identified another effect of thermal fluctuations in a mir- 
ror - due to thermorefractive factor 0i,h = dni^h/dT in 
the layers of the coating, the phase of the reflected way 
should be also proportional to these thermal fluctuations: 
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This value of /3e// is valid only if the topmost covering 
quarter- wavelength layer (cap) has lower refractive index 
ni < Uh and thickness dc = A/(4n;) (A is the wavelength) 
and incorrect in the opposite case, when dc = 0. More 
general case for an arbitrary dc is analyzed in Appendix 
A and discussed below. The existence of excess ther- 
morefractive phase noise was earlier predicted [17] and 
measured [18] in optical fibers. This noise was also mea- 
sured in high-Q optical microspheres [11]. 

VI. Substrate photothermoelastic noise. Noise 
may be produced not only by intrinsic thermal fluctua- 
tions but also by fluctuation of absorbed power, heating 
the mirror [4]: 
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where luq is optical frequency and Sats is spectral density 
of the absorbed power. For the case of shot noise in 
the absorbed power Pabs^ single sided spectral density 
Sabs = 2hujPabs- More general case of arbitrary power 
fluctuations was analyzed in [19]. 



VII, VIII Coating photothermoelastic and pho- 
tothermorefractive noises. Analogously to IV and V 
extrinsic thermal fluctuations produced by absorbed op- 
tical power should produce coating noise due to ther- 
mal expansion and fluctuations of the refractive index 
in the surface layers, with approximate estimates ob- 
tained in [20], which are verified in Appendix B using 
BGV approach, where low frequency approximations for 
G™"j(a;) are also obtained: 
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IX, X, XI Substrate and coating Stefan- 
Boltzmann thermoradiation noises. In all the above 
derivations Stefan-Boltzmann radiation from the sur- 
face was neglected considering that thermal conductivity 
dominates this process. However thermal radiation as 
a dissipative process produces additional fluctuations of 
temperature applied to the mirror. The spectral density 
of these temperature fluctuations obtained by van Vliet 
is given in [15, 16]. As this noise is applied to the surface 
practically in the same way as photothermal noise, we 
can use expressions obtained in VI, VII, VIII to estimate 
thermoradiation elastic noise in the substrate and in the 
coating and thermoradiation refractive noise by simple 
substitutions 
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All the equations for the noises above are given or nor- 
malized for the case of infinite semi-space and adiabatic 
frequency range and show the ceiling of the noises in the 
range ^ >/>■%(£> = n/ipC) is thermal diffusivity). 

This range is mostly valid for the operating frequencies of 
LIGO. Corrections for higher and lower frequencies and 
finite size mirrors where needed, which could be essential 
for other applications with methods of accounting them 
for some cases may be found in [21] and [22]. 

All the noises numbered above are plotted on Figure 
1 for the parameters of LIGO mirrors given in Appendix 
E. It is clearly seen that currently only Brownian and 
coating thermoelastic and thermorefractive noises are es- 
sential. It is appropriate to note however that the Brow- 
nian noises are determined by the loss factor which is 
gradually improving with the development of technol- 
ogy and has not reached as it looks a fundamental shelf. 
This is especially true for the case of coating Brown- 
ian noises where the loss factors in thin films are sev- 
eral order higher than in the bulk. It is ciuitc possible 
that upcoming progress in the technology of dielectric 
film deposition can significantly improve the situation. 
It is worth to note in particularly recent significant im- 
provement of tantala layers losses by addmixing of ti- 
tanium oxide [14]. Another improvement is possible by 




using multilayer coating consisting with unequal layers 
and in particular generalized quarter-length layers where 
high-refractive layer has smaller and low refractive has 
larger optical thickness, preserving half a wavelength in 
the sum [23]. As it is high-reflective material that has 
significantly higher losses and larger thermal expansion 
it becomes possible to lower Brownian surface and ther- 
moelastic noises optimizing relative thickness and com- 
pensating decreased reflectivity by additional layers [24] . 

Thermorefractive and thermoelastic noises unlike 
Brownian noises are determined by basic constants of 
the material which could not be easily modified. These 
noises can become dominating at frequencies higher than 
1 KHz in projected LIGO mirrors. One may notice that 
these two noises having the same order of magnitude and 
the same dependence from frequency and beam waist are 
produced by the same mechanism - temperature fluctu- 
ations as well as other pairs VII- VIII and X-XI where 
these fluctuations of temperature are produced by exter- 
nal source or by surface radiation. Surprisingly a natural 
question about correlation of these noises has not been 
explored so far. 

III. COMPENSATION OF THERMO-OPTICAL 
NOISES 

H.J.Kimble recently proposed a magnificent idea [25], 
that surface fluctuations and refractive index change pro- 
duced by strain could compensate each other in specially 
designed multilayer coating. In principle these strain pro- 
duced fluctuations should be added to the "noisy zoo", 
however strain fluctuations should be significantly sup- 
pressed near the free surface. As Brownian fluctuations 
are currently the main limiting factor in LIGO mirrors, 
the realization of this idea can radically improve the sen- 
sitivity. 

It is shown below that the same idea can be very eas- 
ily applied to the surface noises produced by tempera- 
ture fluctuations, namely to thermoelastic (IV) -I- ther- 



FIG. 2: Schematic of thermal noise compensation 

morefractive (V) [26] as well as pairs of photothermal 
(Vl-fVII) and thermoradiation (X-I-XI) noises. 

Phase fluctuations produced by thermal refraction are 
equivalent to those produced by effective surface fluctua- 
tions Sxeff =54) — XPgffST. I have removed a confusing 
minus in the definition of /3e// a-s compared to BGV [5] 
to have Pef f positive and have the same sign as (3 = 
which is positive in most of optical materials). The ef- 
fect of thermal refraction leads to lengthening of optical 
thickness nd and moves the effective surface from which 
the beam is reflected deeper in the mirror away in the 
direction of incoming beam. At the same time thermal 
expansion moves the surface of the mirror in the oppo- 
site direction toward the beam. This speculation and an 
observation that the value of thermorefractive and ther- 
moelastic coating noises are close enough (Fig.l) gives 
hope that these two effects can compensate each other if 
some parameters of the coating are tweaked. 

The simplest model (Fig. 2) can illustrate the effect of 
thermal noise compensation. Let we have 2A^ layers de- 
posited on unmovable substrate and let the wave is re- 
flected from the outermost layer only. The system of 
equations for the stationary amplitudes of circulating and 
reflected waves 

= -raoe^'''^ + r^aoe^'^'''''^'+'^^ (12) 

I use cx e-«'^*+*'=^ for a plain wave incident on the mirror, 
where k — 2ti/\ is the wave number. With the reflec- 
tion coefficient r — the system has the following 
solution: 

-p _ ar _ 2ikL 6*"*" ~ ^ 
^ — ^ ^ 1 i6' 

ai 1 — re^" 

<j)(. — 2kdc'nc (13) 

Now if due to the change of temperature the index of 
refraction and the length are changed Uc ^ Uc + /3AT, 
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dcTic dcnc{l+(3/ncAT+aAT), and L L — ad^ AT, 
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One sees that though the last term in brackets is os- 
cillating, its average linear dependence on dc allows to 
compensate the first term caused by thermal expansion 
by the thermal dependence of the refraction index with 
appropriate choice of the cover layer thickness. 

In Appendix A it is shown that the same situation is 
also valid for multilayer coating with varying outermost 
low-refractive layer. 



eff,N^oo 



(15) 



1 TTTifiPi + Ph) + A (0c - sm(j)c)inl - nf) 
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The same type of expressions for limited number of layers 
and unequal layers are too bulky for practical purposes 
and it is easier to calculate them numerically by recurrent 
procedures. 

Fig. 3 plots the dependence of absolute value of Peff 
from the thickness of the covering layer. Horizontal 
lines show the equivalent level of thermoelastic param- 
eter 2ac(l + o's)/X. The crossings of the /3-curve with 
horizontal a-lines provide the points of compensation (at 
low frequencies). Solid curves are calculated for the QWL 
coating with 21 layer, while dashed curves stand for op- 
timized unequal 30 double layers with dh = A/8 and 
di — 3A/8, providing approximately the same low trans- 
mittivity ~ 10~^. The dependence of transmittivity (ne- 
glecting optical losses) on the thickness of the topmost 
layer is shown in the bottom of Fig. 3. In most coating 
designs they choose the cap to have zero zero or X/2ni 
(protective cap) thickness giving the best reflectivity, i.e. 
minimal transmittivity. However degradation of reflec- 
tivity can be more than compensated by one additional 



pair of deep layers if the noise issues are considered as 
more significant. It is possible to tweak simultaneously 
both high- and low-reflective top layers. As thermore- 
fractive noise at dc = is somewhat smaller than ther- 
moelastic noise in a wide range of frequencies both for 
equal QWL and for layers with optimized thicknesses, 
noise compensation benefits from optimized levels in the 
same way as thermoelastic noise. 

Noise cancellation can not be absolute. The effective 
volumes of thermorefractive and thermoelastic layers are 
significantly different as the wave is mostly reflected from 
a few outermost layers where the optical power decays 
exponentially and the thermal expansion is provided by 
the whole structure of the coating. The analysis of the 
level of possible compensation is given in Appendix C. 

The results of this analysis are shown on Fig. 4 for the 
case of optimized layers {dh = X/{8nh), di = 3A/(8n/)) 
and cover layer optimized for the frequencies of opera- 
tion 10, 100 and 1000 Hz. For comparison thermorefrac- 
tive noise in QWL layers leveled at thermoelastic noise 
is plotted. One can see that the proposed simple idea of 
suppression of thermorefractive and thermoelastic noises 
allows to lower their combined level up to 3 times in 
the vicinity of lOOHz-lKHz range. If the noise is com- 
pensated for zero frequency, the combined noise becomes 
frequency independent at low frequencies. At high fre- 
quencies, when the length of thermal diffusion becomes 
comparable with the thickness of the coating, the spec- 
tral density of thermoelastic noise falls more rapidly than 
that of thermorefractive noise and compensation though 
at lower level becomes worse. Roughly speaking, simple 
leveling of alpha and beta parameters should work for 
low frequencies below 100 Hz and QWL covering layer is 
close to optimal, while for higher frequencies very thin or 
half-wavelength covering layer, providing minimal /3e// 
is better. It is appropriate to note however, that the 
method of analysis used in Appendix D becomes too ap- 
proximate in this case and more accurate calculations are 
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FIG. 5; Multilayer coating 



required for this case, which are forthcoming. 

It is instructive to look also whether it is possible to 
compensate with thermorefractive effect not only surface 
but substrate thermoelastic noise. The analysis is given 
in Appendix D. It is shown that because of quite differ- 
ent effective working volumes of these effects the level 
of compensation is negligibly small for LIGO mirrors 



0.015/ y/f opt), moreover in LIGO 
the combined even compensated effect in the coating is 
still higher than the effect in the bulk for the operating 
frequencies, however this compensation may be essen- 
tial in more traditional interferometers with significantly 
smaller beam waists. 



Conclusion 

Summing up the results of the presented analysis we 
see that the list of essential noises in high-reflection mir- 
rors by simple tweaking of a multilayer coating may be 
noticeably shortened and we return to the situation of a 
decade ago when only Brownian noise was considered. 
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APPENDIX A. MULTILAYER COATING 

Let we have a multilayer coating, consisting of 2N al- 
ternating layers (or N double layers) on a substrate, with 
the first being a cover layer and the substrate denoted by 
2iV-|- 1 (Fig.5). 



The method of calculation of properties of a multi- 
layer reflection coating is quite simple and computation- 
ally straightforward (see for example [27]), however not 
without pitfalls. Though several analytical expressions 
are obtained below, for different practical applications it 
is usually more appropriate to do numerical anaysis. 

Taking the dependence of forward and backward plane 



waves E,H±(z) as 



E,H{z) = E,H+{z) 



E,H^ (z) one introduces two values - amplitude coeffi- 
cient of reflection r(z) = E-{z) / E+{z) and impedance 
Z{z) — E{z)/H{z), connected by the following equa- 
tions: 



Z{z) 

r(z) 




(16) 



As tangential components of electrical and magnetic 
fields are continuous on the boundary, the same is true 
for the impedance Z. From the other hand, F jumps on 
boundaries but change phase only in homogeneous media 
r(z + d) = F(z)e-*2'="'^ = r{z)e-"l'. In this way, starting 
from the rightmost layer, where in the substrate Fj = 
and Zs = ijs, moving layer by layer from right to left, we 
can calculate Tout- 



F. = 

^out 

F 



F 



27V+1 
^3 + 1 - 



Zj+i + rij 



= 



1 + r^j+iFj 



,i<Pj 



F, 



oc 



F, 



1 + TocTc 
_ V3+1 - Vj 



Vj+i + nj + Uj+i 



(17) 



The most frequently used type of multilayer coating is 
the one consisting of alternating high and low refractive 
quarter wave layers (QWLs) with n2j = Uh, ft 2^+1 = ni, 
dj = (f)j = 2nkdj = tt. This type of coating produces 
the highest reflectivity for the given number of layers. It 
is easy to show that in this case Zj — -'^'-'^ 
different cases 1) when the cover layer n.c = ni is also 
a QWL ~ low refraction output layer, and without the 
cover layer (dc = 0, high refractive output layer) one 
gets: 



1 



eus 



outA 



1 + eus 
1 — eus/n 



1 + ens/nf ' 



(18) 



where e — [rii/nhf'^ . It is essential that when N 
00 and ni/rih < 1 both cases give |Fout| = 1 but with 
different signs. 
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When due to some homogeneous fluctuation in the mir- 
ror coating thicknesses and refractive index change, the 
output phase wifl also change. To calculate this phase 
fluctuation in a general case of arbitrary cover layer we 
use the model of infinite coating [5]. 

Let we have a multilayer coating with infinite number 
of alternating layers of optical thickness nidi and nhdh 
and in some layer nh near its left boundary we have the 
reflectivity Too- Now we add from the left two other 
layers n; and n^^. 



1 + TihToo 1 + ri/jFoo 



hi ^l ^2ikdhnh _ 



1 - TihTi 



e*"^" (19) 



As again we have the infinite number of the same type 
layers then = Too and we can solve the system for the 

F • 



p2 , 

^ on I 



Tih{l - 



We may find the solution of this equation for the gen- 
eralized QWL layers with 2kni^hdi,h = tt ± ipg, e*'^'''* = 

Substituting T^o = e''^^ in (20), we obtain the follow- 
ing equation: 

sin + sin((^oo + V's ) + ^ih sin ipg = 0, (21) 
and after some manipulations we find: 

L = - ^ - arcsm ( Tin sm 



2 
2 



" - ' arcsin [Tih sin ) + tt (22) 



ni ^ ni + Piu 
rih^ nh+ PhU 

i 

7h = - 



(23) 



4^/1 - Tf,^ sin^ '-^nmhTih cos ^ 



niP,{n - Vs) (r^, cos ^ - ^1 - rf, sin^ ^) 



+nhf3i{7: + iPg)il-Tl) 



(24) 



Now we add a cover layer with dc and Uc = ni and 
calculate Tout- 



^out 



ni-nh + im + n/,)rg^ ^2^kd.n,, 
ni+nh + {ni - n/j)r^ 
l-ni + {l+ni)Te 
1 + -h (1 - n;)rc 



(25) 



In particular case of vanilla QWL layers under homo- 
geneous fluctuation of temperature u we have: 



rih 
ni 



(26) 



In principle wc may take into account here also thermal 
expansion as — Pi^h + ni^hO-i^h, however in most of 
the materials used in multilayer coatings a <C p. That is 
why we omit it here retaining thermal expansion only in 
thermoelastic noise where the effect is multiplied on the 
number of layers. 

In the flrst order on u: 



1 oo = -1 + W?j7r^ 



(27) 



(28) 



With the cover layer we obtain 
To = ^^^^^e'^'il + in^Piu) 

= -(l + ^.n./- + ^'[^-f^(^"^/-')])e-^^ 

Tout = C(l+4z7r/3e//«) 
pO _ ■'^f ^ 1 + + 1) COS 4>c + "iini sin cpc 



nf + 1 + {nf — 1) cos (f)c 
_ 1 7rnj^(A + Ph) + A('/'c - sin(/)c)(n^ - nf) 



2-K {nl - nf){nf -|- 1 + (nf - 1) cos^c) 



If rfc = then we have a QWL multilayer coating start- 
ing from high index. 



Lout = -(1 + iT^— 2^) 



'/i ^ 



Peff,H - 7-2 -2 



4 n;; - 



(29) 



If dc = di, then we have QWL multilayer coating start- 
ing from low index. 



nfl3h + nlPi 



1 nfHh 



nt - nr 



(30) 



This expression coincides with the one initially found 
by BGV in [5] and disproves correction proposed in [6]. 
Recently this expression has been also confirmed by Pinto 
[28]. Note that the expression for dc = could not be 
simply obtained by indexes I and h reversal. The pitfall 
here is that these results belong to different branches of 
solution with ^ 1, Z — > 00 for the first low refraction 
layer and T^o ^ — 1, — > for the first high refraction 
layer. 
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From the equation (29) it follows that /?<;// may be 
considered as a sum of fluctuating term and a term lin- 
early changing with dc- Note that for a positive variation 
of temperature and positive phase shift and hence 
/3e// is also positive. 

The equation for a a finite number of 2N layers may be 
obtained in the same way as (18) but keeping first order 
small terms: 

^out,2N = 

{nf - 1)(1 + e^) + {nf + 1)(1 ~ e^) cos 0^ 
{m + e)2 + (1 + me)^ + cos(Ac(n2 _ i)(i _ e^) 

2;)/(l — c') sill (J,. 



+1 



{ni + e)2 + (1 + nisf + cos0c(«? - 1)(1 - e^) 



4:s[Nl32ni/nh - {N - cos(^c)A] 



^.■ J!i^-^) inf - enl) + /3^(1 - efnl 



+2i/J;>e(l - e') - 2iA(l + £') sin</)e] } (31) 

A = (1 - e^) (nf + 1) + cos 0c(nf - 1)(1 + £^) 
-2iesin^c(ni^ - 1) (32) 

This equation transforms into (29) for e — > 0. 



APPENDIX B. PHOTOTHERMAL NOISE IN 
THE COATING 

To calculate photothermal effect in the coating we ap- 
ply to BGV model [4]. The power is absorbed in a thin 
layer dr of the mirror where the internal intensity decays 
e times: For the QWL coating 



dr = 



X{ni + Hh) 
Snirih \n{nh/ni) ' 



(33) 



volume: 



u{lu) 



ndrw'^ 
2w{uj) 

2w{uj) 



-2{v^+Z^)/w^-z/dr. 



u{r,u)) d^r 



/// 



g-«;=(feJ+fe^)/4 



{Dk'^ + iijo){l - ik^dr) (27r)3 



pC J J Dkl 

O-oo 
W{U!) 



, ,.,,2 



iuj (27r)2 
VTrift^e*™ [1 - erf (v/i6^)], 



2D 



(36) 



Finally as the temperature fluctuations produce phase 
fluctuation equivalent to x{uj) = (3eff\u{uj) we obtain 
the spectral density: 



Qcoat 



•eff) 



5 
462" 



sur. 



1 



bn, > 1 



(37) 



For very high frequencies when br = ^jjf ^ 1 this ap- 
proximation breaks down and should transform into the 
substrate photothermal noise with substrate parameters 
substituted by parameters of the coating. 

This expression coincides after substituting Sabs with 
that obtained in [20] for 3> 1. For practical purposes 
it is possible to match a simple numerical approximation 
for J which has maximum error less than 5% for all 



In case of the silica-tantala QWL coating and A = 1.064 
mkm dr — 0.43 mkm. As dr <^w,rT = al\f^ the model 
is the following (A.l): 



du 



w{t) 2 



at pL 'KW^ 



(34) 



where w{t) is the fluctuating part of the absorbed power 
w{t) = Wabsit) — {Wabs) which for the shot noise has the 
spectral density Sabs = ^^oi^abs)- The solution is: 



+ 00 



"(i-''^) = —FT- / / / KT^rr—. — (35) 



pC JJJ Dk^+iLo (27r)3 

— oo 

Temperature fluctuations averaged over the effective 



wbu, - 2.20bT + 8.9862 
1 + TAAbyj - 3.456^/^ + 8.9862 



(38) 



To calculate the thermoelastic effect of absorbed power 
in the layer djv one should at first solve for thermoelastic 
problem with the given fluctuating temperature to find 
surface fluctuations which should be averaged over the 
beam spot. This problem has been solved in [6] (A. 15) for 
the Langevin fluctuations: 



x{uj) = 
2adN{l + (j) jjj< 



(39) 



+ 00 



u(r, w) d^r 



Comparing this expression with (36) we sec that bound- 
ary conditions on free surface produce more than doubled 
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effect as compared to naive x{uj) = adNu{uj). Hence we 
obtain 



surf ' 



(40) 



This estimate is two times smaller (for G"""^^ = 1) than 
that obtained in [20] where also Poisson correction is ab- 
sent, however one should note that in real coating with 
very inhomogeneous initial temperature distribution over 
3-4 layers and inhomogeneous a and pC , coefficients may 
sufficiently differ from these rough approximations and 
from averaging adapted from surface thermoclastic cal- 
culations in [7] high frequency limit of this approxima- 
tion is reached at significantly lower frequencies, when 

bd = ^^5^ 1, where the effect is transformed into the 
bulk photothermal effect of BGV [4] in the effective layer 
material. 



APPENDIX C. THERMOREFRACTIVE AND 
THERMOELASTIC NOISES 

The calculation of cross-correlation effects could not 
be easily obtained using simple fluctuation-dissipation- 
theorem approach. In this case the Langevin fluctuating 
forces are valuable. Using the previous BGV results the 
spectrum of the fluctuations of temperature in the mirror 
may be found as: 

-t-oo 

.(".r)=///S^^ (41) 



Dk'^ + iLU (27r)3 



{F(k,Lu)F*{k',u;')) = 
(27r)4^^fc2^(k-k')5(a;-a;') 

Thermorefractive fluctuations in a very thin effective 

layer where the reflection mostly occurs are propor- 
tional to the exponent- weighted averaging over the effec- 
tive reflective layers dr- 



(42) 



+ 00 +00 



rr r 

ndrw'^ J J J 



u{uj, r) 



-oo 



-/3e//A 



{D,k^ +iLo){l-ik^dr) (27r)3' 



where fc^ = ky + k'^. As the layer is thin the main con- 
tribution to the fluctuations of temperature in the layer 
is produced by substrate. As it is shown in [7] the con- 
tribution of the layer itself is ^ iodj^ / Dc times smaller. 
That is why all the parameters under integral are related 
to the substrate. The effect of the layer becomes ap- 
parent due to thermal expansion and fluctuation of the 
refraction index. 



The calculation of the effect produced by thermoclas- 
tic effect in the whole multilayer coating on the substrate 
is a more difficult problem. It requires the solution of the 
thermal and elasticity problems with a subsequent aver- 
aging of the surface displacement. However as the spec- 
tral density has been already accurately calculated [6, 7], 
one may in the flrst order on rfjv use analogous averag- 
ing with effective thickness d* ~ djv fitted to the known 
results. An accurate straightforward calculation for inho- 
mogeneously distributed noise sources may be performed 
using the approach of Green's functions for the Langevin 
sources developed by van Vliet [15, 16]. 

\t) ^ (43) 



Y'coat ( 



+ 00 



(277)3 



-fe1 



d^k 



{D^k"^ + iLo){l - ika^d*) 



We obtain: 



Qcoat Qcoat 

•^TR ~ '->TER — 



SDskeT^ fff 2a,dN{l+(Ts) 

k^e 



(44) 



1 — ikr.d* 



1 %J^rndr 



— OO 



Qcoat 



1 + fc2d*2 £>2fc4 + ^2 

— CSO 

^ 8kBT^aldl{l + a,)2 1 - ^6 + 6^ 

~ 7r3/2,«2y^;^/l/2 1 + ^,4 

b=^/cod*'^/Ds 



, (45) 



Prom (6) it follows that \/26 ~ ^w^^ ^^^^^ ^* - 



dN 



3R-'-l 
3R • 



Qcoai 
•^TR 



2kBT^0iff\ l-V2c + c2 

3/2^2^^/1/2 1 + ' 



(46) 



c= ^/ujdl/D 



QCoat 
^TER 



SkBT^acdNjl + cTs)peffX 1 - V2(6 - c) + b{b - c) 

We may now choose a frequency of interest fopt and 
optimize for it the combined thermoelastic and thermore- 
fractive noise. Considering c<^b: 

f3eff,opt)^ — 2acdN{^ + Crs)B{fopt) 



B{f) 



l-V2b + b'^ 

1 + 64 



(47) 
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In the calculations I assumed that in the area of in- 
terest dr <C <C w and k± <^ k^. The approximation 
can possibly break down when ^ ^ 1 which is attained in 
silica-tantala coating with 20 double layers at / ~ IkHz. 
In this case thermal sources in the coating itself should 
be additionally considered. 



Now optimizing /3e/ / at a frequency of interest we find: 



S^,min = S^t[l-2j-^] (49) 



APPENDIX D. SUBSTRATE AND COATING 
THERMOELASTIC NOISE CORRELATION 



The effects of thcrmoelastic substrate and coating ef- 
fects arc the following: 



2a,(l + £7,) 



_fk_ f(k,C^) k^ i^t-klw'/8 

(27r)3 D,/e2 + fc2 



(277)3 D.k^ + iu 

— OO 

SDskeT^ 



2as{l + as)^-PeffX 



-d^k 



Qcoat I Qsub Q 

, ^ 4kBT^al{l + asrn, 
2kBT^I3lff\^ 



QCoat 



Ster = 



■ ff r-''" rill kjdk^dk^ 

J J D2fc4 + W2 

— CO 

^ 8kBT^^as{l + cr,)l3effX 



TT' 



APPENDIX E. PARAMETERS USED IN 
NUMERICAL ESTIMATES [7] 



General parameters: 



T = 300K A = 1.064 X 10-6 m 

Wo = 1-77 X lO^s s-i w = 0.06m 
Fused silica substrate parameters: 



as = 5.1 X 10"^ K-\ = 1.38J/(ms K), 

p = 2200 kg/m^ C = 746 J/ (kg K), 

Y, = 7.2 X 10l°N/m^ a = 0.17, 
d> = 5x 10-9 



Thin film parameters (Silica and tantala): 



ni = 1.45 
ai = Us 

Kl = Ks 

Hh = 2.065 

ah = 3.6 X 10-*^ K-\ 
Kh = 33 J/(m s K), 



4>i = 1.0 X 10-4 
f3i = 1.5 X lO-^K-^ 

Pi = Ps, 

<i)h = 3.8 X 10-4 
A = 6 X 10-5 K-\ 
p,, = 6850kg/m^ 
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